Spectral Properties of 5-Plutonium: Sensitivity to 5/ Occupancy 
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By combining the local density approximation (LDA) with dynamical mean field theory (DMFT), 
we report a systematic analysis of the spectral properties of 5-plutonium with varying 5/ occupancy. 
The LDA Hamiltonian is extracted from a tight-binding (TB) fit to full-potential linearized aug- 
mented plane-wave (FP-LAPW) calculations. The DMFT equations are solved by the exact quan- 
tum Monte Carlo (QMC) method and by the Hubbard-I approximation. We demonstrate strong 
sensitivity of the spectral properties to the 5/ occupancy, which suggests using this occupancy as 
a fitting parameter in addition to the Hubbard U. By comparing with photoemission data, we 
conclude that the "open shell" 5/ 5 configuration gives the best agreement, resolving the controversy 
over 5/ "open shell" versus "close shell" atomic configurations in 5-Pu. 

PACS numbers: 71.20.Gj, 71.27. +a, 75.20.Hr, 79.60.-i 



I. INTRODUCTION 

From a consideration of its condensed matter physics 
properties, crystal structure, and metallurgy, plutonium 
is probably the most complicated element in the peri- 
odic table jii2iMi5i£ including a phase diagram with six 
allotropic phases. The low-temperature monoclinic a- 
phase is stable up to 395 K while the face-center-cubic 
(fee) (5-phase is stable between 592 and 724 K. Further- 
more, in stabilized alloys the a to <5 phase transforma- 
tion of Pu has a significant volume expansion, with the 
(5-phase 25% larger in volume than the a-phase. This 
behavior is related to the special position of Pu in the 
periodic table, which is at the boundary between the 
light actinides that have itinerant 5/ electrons and the 
heavy actinides with localized 5/ electrons. In this situ- 
ation the electrons are in a very strongly correlated state 
where even the best conventional LDA band-structure 
calculations cannot predict its unique properties, e.g., the 
<5-phasc volume^ in the experimentally observed non- 
magnetic state. This failure has stimulated numerous at- 
tempts to include additonal electronic correlation effects. 
Although LDA theory has repeatedly claimed 9 - ' 10 : 11 ' 12 : 13 
that the thermal expansion in Pu is a consequence of 
magnetism, this is in striking contradiction with exper- 
imental data.— Several research group o 15 ' 16 ' 17 have ap- 
plied the LDA+U method to include more /-/ correla- 
tion energy. By adjusting the on-site Coulomb repulsion 
energy (U) appropriately, it has been possible to fit to the 
experimental (5-phase volume. The LDA+C/ calculations 
also indicated an instability of 6-Pu toward an antifcrro- 
magnetic ground state. More recent calculations based 
on the "around-mean-field" LSDA+l7 including the spin- 
orbit interaction^ or the fully-localized-limit LSDA+£7 
quenched the spin polarization through an w / 6 configu- 
ration, and obtained the correct (5-phase volume due to a 
weak exchange interaction.— To predict the correct vol- 
ume, a mixed-level model has also been proposed by one 
of us together with co-workers^ In this model, four of 



the five 5/ electrons in the (5-phase were constrained to be 
localized and were not allowed to hop from site to site, 
but could hybridize with the conduction electrons. A 
major advance came when a new calculation scheme was 
proposed to merge LDA-based methods with DMFT . 21 : 22 
DMF T 23 ' 24 is a many-body technique that is able to 
treat the band- and atomic-like aspect simultaneously 
when applied to plutonium. Within the LDA+DMFT 
approach ) 22 : 25 the origin of substantial volume expan- 
sion was explained in terms of the competition between 
Coulomb repulsion and kinetic energy. Thus, two dra- 
matically different pictures for the non-magnetism of <5- 
Pu have emerged. In the "5/ 6 " picture ) 18 : 19 ' 26 : 27 by in- 
cluding the spin-orbit interaction, one starts with a closed 
/ 5 / 2 atomic subshell fully filled with six electrons while 
the / 7 / 2 subshell is empty, making Pu magnetically in- 
ert. In the "5/ 5 " picture ; 25 : 28 one starts with an open 
/ 5 / 2 atomic subshell filled with five electrons, resulting 
in a magnetic moment that should be screened by the sd 
valence electrons. Therefore, more studies are necessary 
to determine which picture will prevail. 

In this paper, we report a systematic LDA+DMFT 
study of the spectral properties of <5-Pu with varying 
5/ occupancy. Throughout this work, the LDA part of 
Hamiltonian is determined from a new TB fit to FP- 
LAPW calculations. The DMFT equations are solved 
using quantum Monte Carlo (QMC) simulations as well 
as the more approximate Hubbard-I method. These pro- 
vide an accurate characterization of spectral properties of 
<5-Pu. It is found that the 5/ spectral density of Pu is very 
sensitive to its occupancy. Good agreement is found with 
photoemission spectroscopy (PES) measurements when 
about five electrons occupy the / 5 / 2 subshell, in support 
of the second picture. Other recent work appears also 
to be in support of this picture. Measurements of the 
branching ratio in bd to 5/ transitions favor a 5/ count 
closer to 5 than 6r£>2P. while subsequent LDA+DMFT 
calculations^ using a vertex-corrected one-crossing ap- 
proximation for the auxilliary impurity problem yield re- 
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suits for the branching ratio in agreement with experi- 
ment^ and suggest a 5/ occupancy of 5.2. Similarly, spe- 
cific heat calculations carried out using an LDA+DMFT 
method with a perturbative T-matrix and fluctuating ex- 
change approach to the auxiliary impurity problem also 
favor the ~ 5/ 5 configuration.— 

The remainder of this paper is organized as follows: 
Theoretical methods are given in Sec. [Til numerical re- 
sults for the 6-Pu are presented in Sec. IIII1 and then a 
summary follows in Sec. IIVI 

II. THEORETICAL METHODS 

Most LDA+DMFT methods to date have been imple- 
mented in a basis of linear muffin-tin orbitals (LMTO) 
within the atomic sphere approximation (ASA), or in its 
full potential version. Here we use an alternative ver- 
sion of the LDA+DMFT method, which is based on a 
recently developed TB theory^ In this representation, 
the full second-quantized Hamiltonian is written as: 

+ -^- ^ l ifjmjfl'ifi'm' j ■ (1) 

Here k are the Brillouin-zone wavevectors; i are lattice 
site indices; I is the orbital angular momentum; j is the 
total angular momentum; rrij = —j, —j + 1, j — 
n l f ]mj = clj jm .Cifj mj . The relevant orbitals for Pu are 

7s, 6p, 6d, and 5/, and so the matrices H°(k) are 32 x 32. 
They are given by 

H°(k) = ff LDA (k) + ( e/ -4 DA )J / . (2) 

Here the matrix If is zero except for l's along the 14 /-/ 
diagonals, and 

£ / DA = I^vE Tr ^ LDA ( k ) f /]= (3) 

k 

where N is the number of k points in the Brillouin 
zone. The matrices _ff LDA (k) are orthogonalizcd, single- 
electron Hamiltonian matrices obtained from TB fits to 
the FP-LAPW calculations^ In this study, strict orthog- 
onality is maintained between the TB orbitals (hence no 
overlap matrix need be included). All matrices are cal- 
culated at the experimental i5-Pu volume. 

The Slater-Kostcr tables for the sp 3 d 5 matrix ele- 
ments can be found in standard references ; 33 ^ 4 and we 
have used an extended formalism for a unified treat- 
ment including additional matrix elements involving /- 
electrons.— Typical TB calculations are then reduced to 
using TB as an interpolation scheme; the matrix elements 
are determined by fitting to ab-initio calculated quanti- 
ties such as the total energy and band energies. In this 



study TB inter-site parameter values are evaluated at 
inter- atomic distances out to the fifth nearest neighbor, 
resulting in 100 inter-site parameters and 4 on-site pa- 
rameters. The effect of the spin-orbit interaction is in- 
cluded as a perturbation, with its parameters kept fixed 
across the bands, although evaluated at the most impor- 
tant energies, namely at the respective centers of gravity 
of the occupied state density for each orbital type^- The 
spin-orbit coupling adds 3 more parameters to the TB fit, 
resulting in a total of 107 parameters (parameter values 
are available on request from the authors). The quality of 
the present TB fit for <5-Pu is comparable to that shown 
elsewhere for fee Ui22 

As customary in LDA+DMFT calculations, Eqs. Q} 
12]) presume that the LDA results are sufficient for both 
the off-diagonal elements of ff°(k) (hybridization) and 
the diagonal elements of the assumed uncorrelated spd 
electrons. For the correlated / electrons one needs of 
course the Coulomb interaction as seen in Eq. (UJ, but 
in addition, since the LDA site energy defined by Eq. © 
includes Coulomb contributions, it must be replaced by 
the true or bare site energy Ef as in Eq. (J5J). Both Ef and 
Uf may be estimated from the dependence of the LDA 
site energy £y DA (n,/) on / occupation rif via LDA con- 
strained occupation calculations , 37 i 38 and comparison to 
spectroscopic data for ions suggests uncertainties of ±1 
eV in such results for both quantities^ Typical values 
of Uf are around 4 eV for <5-Pu, as we also find, and this 
value is assumed throughout the present paper. Given 
the sensitivity of n/ to £/, on the other hand, we have 
chosen to adjust the parameter ej in the present work 
as a natural and convenient means of exploring the sen- 
sitivity of the (5-Pu spectra to 5/ occupancy^ This ad- 
justment amounts to varying the double counting term 
in conventional LDA+C/ or LDA+DMFT terminology. 

Within the DMFT, the lattice problem Eq. (TJ is 
mapped onto a multi-orbital quantum single impurity 
problem subject to a self-consistency condition: 

g-\iu) n ) = G^(iw„) + t(iw n ) . (4) 

Here Q(ioj n ) is the Weiss function, S(icJ n ) is a k- 
independent self-energy, and the local Green's function 
is defined as G\ oc {ioj n ) = G-k(iu) n )/N, where the lat- 
tice Green's function reads 

G k (w„) = + fi)I - H°(k) - tiiun)]- 1 , (5) 

with / the 32 x 32 unit matrix and the chemical poten- 
tial. Since we have restricted strong correlation to the / 
orbitals only, it is customary to take S to have nonzero 
elements only within the 14 x 14 /-/ block, which per- 
mits a similar reduction in the quantum single-impurity 
problem itself. Furthermore, we neglect the off-diagonal 
elements in the self-energy and, by ignoring the crystal- 
field effect in the impurity problem, treat the 5fj =5 / 2 
levels as one 6- fold degenerate level and the 5/j = 7/ 2 as 
another 8-fold degenerate level. The only nonzero /-/ 



block self-energy reduces to 



,{iu) n ) = Y,y(iw n )Sjj>5, 



(6) 



and we need only find functions for the two spin-orbit 
states j = 5/2 and j = 7/2. Note that when no spin- 



orbit coupling is included, Sgy 2 (iti; n ) 
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respondingly, the local Green's functions is dictated by: 
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(7) 

It should also be noted that we neglect the crystal-field 
effects only at the auxiliary impurity problem, all such 
effects are still retained in H°(k), where hybridization 
dominates anyway. 

To solve the auxiliary impurity problem, we have em- 
ployed both the Hubbard-I approximation and an im- 
plementation of the QMC method that has been used 
to study properties of the compressed lanthanidesi 40 ' 41 i 42 
The former method is more approximate but much faster 
while the latter is more rigorous but computationally 
more expensive. The Hubbard-I approximation is also 
useful in various ways: it is more applicable in the 
large-volume atomic limit, where it gives insight into the 
atomic aspects of Pu; it provides a good initial guess 
for the DMFT(QMC) simulations; and it makes it pos- 
sible to examine the full /-/ Coulomb interaction (all 
four Slater integrals F k ) and its consequent term struc- 
ture in DMFT(HI) in contrast to DMFT(QMC), where 
it is still difficult to go beyond just the standard Hub- 
bard Uf = F°. Unless specifically noted otherwise, all 
DMFT(HI) and DMFT(QMC) results reported in this 
paper include only Uf in addition to the spin-orbit in- 
fraction. For the DMFT(QMC) calculations, we used 
6000 sweeps per QMC iteration for the imaginary time 
segments L = 112 at a temperature T = 632 K, which is 
about as low a temperature as is practical at the present 
time with QMC, and performed at least 100 iterations. 
To improve the input data needed for the maximum en- 
tropy (MaxEnt) analytical continuation^ at the last it- 
eration, we dramatically increased the number of sweeps 
such that the number of bins, each containing 100 mea- 
surements, is larger than 2L. 



III. RESULTS 

In Fig. [1] we show the spectral density of <5-Pu at 
T = 632 K for various values of the 5/ occupancy n f . The 
left three panels display underlying atomic properties ob- 
tained from the DMFT(HI) calculations, consistent with 
a ground state a\f^ 2 ) + P\f^/ 2 ), where a| 2 +|/3| 2 = 1 and 
5+|/3| 2 . For nf = 6 in Fig. [TJa), one gets a single 



71/ 

lower (/ 5 6 /2 ^/ 5 5 /2 ) and upper (/| /2 -> f! /2 fj /2 ) Hubbard 
band. For n/ = 5 in Fig.[T^c), however, one may still add 
an electron to the j = 5/2 subshell, so that in addition to 
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FIG. 1: (Color) The 5/ spectral density calculated from the 
DMFT(HI) [left panels] and DMFT(QMC) [right panels] for 
8-Pn for varying 5/ site energy at T = 632 K. In all pan- 
els, the spectral density contributed from the j = 5/2 and 
j = 7/2 sub-shells are represented by the black and red lines, 
respectively. The total 5/ spectral density is represented by 
the blue line. The energy E = denotes the position of the 
Fermi level. 
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one has two, spin-orbit split upper 

is/2 ft/2 and fl/2 -> /sVr/a)- Fot 
the mixed valent, nonintegral n/ in Fig.[ljb), the system 
behaves as an ensemble of \a\ 2 fh 2 and |/3| 2 f^ 2 configu- 
rations. In the absence of hybridization, this leads to four 
Hubbard bands in relative position 4J7/, 5C//, 5C// + A/, 
and 6Uf + Af, where Af is the spin-orbit splitting, with 
areas of 5|cv| 2 , 6\(3\ 2 + \a\ 2 , 8\a\ 2 , and 8|/?| 2 , respectively. 
The Fermi level (energy zero) should lie within the sec- 
ond Hubbard band at 5Uf splitting the occupied (6|/3| 2 ) 
and empty (|a| 2 ) parts, where the latter is the small tail 
of the j = 5/2 spectra (black) extending above the Fermi 
level in Fig.QJb). 

For the more rigorous DMFT(QMC) results one sees, 
in the three right-hand panels of Fig. [1] the expected 
transfer of spectral weight away from the Hubbard bands 
into the quasiparticle peak at the Fermi level, which 
DMFT(HI) is incapable of describing. However, in ad- 
dition there also appears to be a shift of the outlying 
DMFT(HI) Hubbard structure closer towards the Fermi 
level, while at the same time the one Hubbard band al- 
ready overlapping the Fermi level [second from left in 
Fig. [ljb)] spreads away from the Fermi level. The over- 
all affect is to give a more smooth and systematic evolu- 
tion of the DMFT(QMC) spectra with growing n/ than 
seen in DMFT(HI). This is confirmed by an indepen- 
dent evaluation of the DMFT(QMC) state density at 
the Fermi level via Aj(0) = (2j + l)((3/n)G j (T = (3/2)M 
^5/2(0) smoothly decreases while A 7 / 2 (0) (and the to- 
tal) smoothly increase for increasing rif over the range 
5 < nf < 6, consistent with the Fermi level passing more 
into the j = 7/2 part of a spin-orbit split peak. Indeed, 
the MaxEnt results in Figs. [IJd-f) exhibit a splitting in 
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FIG. 2: (Color) (a) Experimental and theoretical photoe- 
mission spectrum, (b) The 5/ spectral density calculated 
from the DMFT(HI) with only F° = U f (black line) and full 
Coulomb (red line) interactions for rif = 5.0. (c) The 5/ spec- 
tral density calculated from the DMFT(HI) with full Coulomb 
interactions for nj — 5.2 and scaled/shifted ones for rif — 5.0 
and 6.0. Here E\ = E + 0.87 eV and E 2 = E - 2.29 eV. 

the quasiparticle peak. On the one hand, we suggest that 
it is due to the spin-orbit interaction and the associated 
j-dependence induced into the self-energy. Such splitting 
appears in the spectra for each j as the two channels 
couple via H°(k). A reanalysis using MaxEnt of ear- 
lier DMFT(QMC) calculations^ for Ce and Pr exhibits 
similar behavior, and since this is then independent of 
filling, hybridization dips or gaps can be ruled out. Such 
spin-orbit induced splitting has been experimentally ob- 
served in photocmission experiments for Ce^ 5 - although 
the Hund's rules exchange omitted in Fig. [1] (and dis- 
cussed shortly) will have larger impact for mutli-/ elec- 
tron atoms. On the other hand, we do not exclude the 
possibility that the splitting may be simply an artifact of 
the analytical continuation technique used to obtain the 
spectral function at the real energy axis. Note that we 
find the quasiparticle peak and its associated fine struc- 
ture to both disappear for S-Pu at higher temperature, 
1580 K, consistent with the 800 K Kondo energy that has 
been suggested^ 

A comparison of the DMFT(QMC) spectra with ex- 
perimental photocmission dat a 20 ' 46 is shown in Fig.[^a). 
The spd valence contribution was added to the calculated 



5/ spectra from Figs.QJd-f), and the total was broadened 
to reflect the 60 meV instrument resolution and 15 K 
measurement temperature, as well as the Lorentzian for 
the photohole lifetime determining the natural linewidth, 
including quadratic scaling with binding energy. For the 
experimental comparison, data from the 40.8 eV HcII- 
alpha line were selected to obtain a photon energy range 
where orbitals of interest have similar cross sections^ 
While additional work is certainly needed for an optimal 
comparison with the data, Fig. [Ua) does reproduce the 
trend seen in the DMFT(QMC) results of Figs. Hfd-f), 
where the lower Hubbard band is seen to move to more 
negative energies with increasing rif, which lends exper- 
imental support for a value of rif closer to 5. 

The implications of the comparison in Fig. HJa) do not 
appear to be compromised by omitting Hund's rules ex- 
change and the associated term structure in the theo- 
retical results, where only the monopolc Slater integral 
F° = Uf has been incorporated. This can be tested in 
DMFT(HI), where the black curve in Fig.^b) is the to- 
tal nf = 5.0 5/ spectra of Fig. [ljc), while the red curve 
now treats the full Coulomb interaction taking in addi- 
tion reasonable experimental values for F 2 , F 4 , and F 6 . 
In the electron removal spectra one now sees the stan- 
dard terms (labelled) of the / 4 final state. These agree 
in relative position and intensity with earlier theoretical 
analysis^ and serve only to broaden the lower Hubbard 
band to the negative-energy side, leaving the positive- 
energy side of importance to the comparison in Fig. [2fa) 
relatively unaffected. While this DMFT(HI) test can not 
probe the impact of Hund's rules exchange on the quasi- 
particle peak near the Fermi level, one might anticipate 
a reduced effect after such fine structure has been ap- 
propriately broadened for comparison to photocmission 
data. 

There has been much discussion of a three-peak struc- 
ture within about 1 eV of the Fermi level in the photoc- 
mission of Pu systems: in thin layers of PuSe^ 8 - thin Pu 
layers on Mg 49 ' 50 single crystal PuTe^i and suggested 
in thin film of PuSii.7^ and PuNj2£ The presence of the 
structure in thin Pu layers (a few ML) was interpreted 
as a result of localization effects due to low dimensional- 
ity overcoming the itinerant character of 5f electrons^ 
More recently it has been suggested that such structure 
may be evident in 8-Pu metal itself, 28 ' 52 Such structure 
may also be seen in DMFT(HI) calculations including 
the full Coulomb interaction. In Fig. |2fc), the black 
curve corresponds to a DMFT(HI) calculation with the 
5/ site energy adjusted so n/ = 5.2. For comparison 
the red (blue/shaded) curve shows the rif = 5.0 (6.0) 
result scaled by \a\ 2 = 0.8 (|/3| 2 = 0.2) and shifted so 
its lowest unoccupied (highest occupied) state is at the 
Fermi level. The three peaks discussed are shown la- 
belled: 6 i?5/2 essentially at the Fermi level, and then 
6 #7/2 and e F 5 / 2 moving below. Given the agreement be- 
tween the DMFT(HI) spectra A 5 . 2 for 71/ = 5.2 with the 
composite O.8A5.0 + 0.2^6.0 (addition not shown), it is 
evident that especially the left two peaks arise from / 5 
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term structure in the final state of 5/ 6 — > 5/ 5 electron 
removal, and that the strength of these peaks might pro- 
vide a direct measure of the 5/ 6 admixture (3 into the 
<5-Pu ground state, and therefore of rif itself. 

There is little doubt that such an analysis is reasonable 
for some Pu compounds, where precise atomic multiple! 
calculations were found in agreement with photoemis- 
sion data in both line positions and intensities^ Such 
agreement includes cases (cubic PuTe and PuSe) where 
the three-peak structure is experimentally observed be- 
low the Fermi level approximately 100 meV. The situ- 
ation provided by a metallic environment as in 5-Pu is 
substantially different, however, and so the localized mul- 
tiplet interpretation is not at all clear. An obvious con- 
cern about the results in Fig. Etc), for example, is that 
the self-energy used in DMFT(HI) is too atomic- like; the 
differences between the DMFT(HI) and DMFT(QMC) 
curves in Fig. Q] add to this concern. 

The nature of the peak right at the Fermi level in 5- 
Pu is clearer. The low-temperature electronic specific 
heat in bulk 6-Pu metal is 64 ± 3 mJ K^moU 1 ^ This 
large, heavy-fermion-like value indicates significant 5/ 
state density at the Fermi level which is consistent with 
the kind of Kondo-likc physics that DMFT(QMC) can 
yield, albeit possibly at temperatures below the present 
work. Corroborating evidence comes from the photon 
energy dependent and temperature dependent photoc- 
mission experiments performed on clean <5-Pu metal sur- 
faces,— which strongly suggest that the first peak 
in the <5-Pu arises from 5/ electrons that are hybridized 
with the conduction electrons, which is again a typical 
quasiparticle feature usually found for materials with an 
enhanced electron mass. 



IV. SUMMARY 

We report a systematic LDA-DMFT study of the spec- 
tral properties of <5-Pu with varying 5/ occupancy. The 



strong sensitivity of the spectral properties to the 5/ oc- 
cupancy and the inherent ambiguity in defining the as- 
sociated 5/ orbitals suggests the need to use 5/ occu- 
pancy as a fitting parameter in addition to the Hubbard 
U for DMFT theories of Pu and other materials. By 
comparing with photoemission data, we conclude that 
an "open shell" 5/ 5 configuration gives better agreement 
with experiment than the hypothesized "closed shell" 5/ 6 
c&sel&l&2&ZL The DMFT(HI) approach, which is ac- 
curate in the localized, atomic limit, gives significantly 
different results from the more rigorous DMFT(QMC) 
method, which confirms the presence of a significant itin- 
erant character in <5-Pu, especially at the Fermi energy. 
The "3-peak" structure seen within about 1 eV of the 
Fermi level in the photoemission spectra of some Pu com- 
pounds as well as metallic thin films, but which remains 
controversial for the bulk S-Pu metal, may provide a mea- 
sure of the admixture of 5/ 6 character in the ground 
state of these materials, and therefore an additional in- 
dication of rif. However, it is not clear from the present 
work whether this structure should wash out in a corre- 
lated calculation for <5-Pu, which rigorously treats both 
extended and intraatomic exchange effects. 
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